ADP-02 80 T519 



On the chirality of quark modes 

Alex C. Kalloniatis* 
Special Research Centre for the Subatomic Structure of Matter, 
University of Adelaide, South Australia 5005, Australia 

Sergei N. Nedelko tt 
Bogoliubov Laboratory of Theoretical Physics, JINR, 141980 Dubna, Russia 
(Dated: February 1, 2008) 

A model for the QCD vacuum based on a domainlike structured background gluon field with 
definite duality attributed to the domains has been shown elsewhere to give confinement of static 
quarks, a reasonable value for the topological susceptibility and indications that chiral symmetry 
is spontaneously broken. In this paper we study in detail the eigenvalue problem for the Dirac 
", operator in such a gluon mean field. A study of the local chirality parameter shows that the lowest 

q ■ nonzero eigenmodes possess a definite mean chirality correlated with the duality of a given domain. 

' A probability distribution of the local chirality qualitatively reproduces histograms seen in lattice 

, simulations. 
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I. INTRODUCTION 



~H ' In a previous paper jlj] we formulated a model which characterises the QCD vacuum by a "lumpy" distribution of 
field strength and topological charge density. For lack of a better name we shall refer to the model as "the domain 
model" . The formulation is given concretely in terms of a partition function which describes a statistical ensemble of 
domains, each of which is characterised by a set of internal parameters associated with the mean background gluon 
field, and the internal dynamics are represented by fluctuation fields. Correlation functions in this model can be 
calculated by taking the mean field into account explicitly and decomposing over the fluctuations. We briefly review 
£NJ . the details and assumptions behind the model in the next section but state here unambiguously that the "domains" 
in question are assumed to be purely quantum in nature. They are not semi-classical solutions of Yang-Mills theory 
, and are not argued to exist as topologically stable classical configurations, rather they seek to characterise the average 
bulk properties of the ensemble of fields that determine the gluonic vacuum. In particular, it is not assumed that 
the topological charge associated with a domain should be an integer. The rationale of such an extremely simplified 
^■i' construction can be understood as an attempt to implicitly incorporate effects of the presence of singular pure gauge 
configurations in the QCD Euclidean functional integral into a practical calculational scheme with a mean field 
• • . description of the QCD ground state. A self-consistent mean field approach requires nonperturbative calculation of 
the free energy as a functional of the mean field whose minima should determine its form, but this is beyond the reach 
of analytical methods. Nevertheless there is an accumulation of semi-qualitative arguments [Q in favour of the ansatz 
5^ ' for the mean field we have chosen. 

In the gluonic sector the model depends on two parameters, a mean field strength per domain B and a mean size for 
domains R, which is sufficient for an adequate description of the pure glue characteristics of the QCD vacuum - the 
gluon condensate, topological susceptibility, and string tension. The Wilson loop in such a gluonic background was 
found to exhibit an area law dependence for large loops. Thus a confinement of static fundamental charges is captured 
by the model; some dynamical gluon degrees of freedom turn out also to be non-propagating. The absolute value of 
underlying average topological charge per domain was determined to be approximately q — 0.15 and the density of 
domains to be as high as 42 fm~ 4 . Although tentative signals of spontaneous chiral symmetry breaking were also 
obtained in Q , a more rigorous consideration of the fermionic spectrum and eigenmodes as well as calculation of the 
quark determinant is required, which was missed in JO . 

In this paper we solve the eigenvalue problem for the Dirac operator for the gluonic background and boundary 
conditions adopted in the model and examine the chirality properties of the eigenmodes. This is a necessary step 
for checking the status of chiral symmetry breaking in the domain model. But in view of recent lattice results this 
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problem is valuable also in its own right. 

There are strong hints in lattice Monte-Carlo simulations at intermediate-range structures in individual gluon 
configurations once fluctuations are filtered out by some means. For example, cooling or relaxation algorithms are 
well established now ||] , and can reveal instantonic like structures after several sweeps of a given lattice configuration. 
However, as these algorithms are designed precisely to locally minimise the action, it is natural they should bring 
objects with integer topological charge into relief. Alternately and more relevant to the present work, low-lying and 
zero-modes of the massless Dirac operator can be used as a probe of long-range gluonic structures || although only 
recently did this become more reliable with lattice fermions with good chirality properties. For example, exact index 
theorems are found to be satisfied on the lattice Q and the zero modes are seen to precisely correlate with instantonic 
structures in the raw lattice configuration, in the absence of cooling However, the exact zero modes of any finite 
volume simulation cannot be those relevant to spontaneous chiral symmetry breaking, rather the discrete spectrum 
of low-lying non-zero modes should, in the infinite volume limit, go over to a continuous band at zero, saturating the 
Banks-Casher relationship Q. Such modes are sometimes called "pseudo zero modes". Low-lying non-zero modes 
with strong signs of chirality in regions of high action and topological charge densities would be a tool for identification 
of the properties of gluonic configurations relevant to chiral symmetry breaking. Indeed, after an initial negative result 
, recent results have emerged showing precisely this: low- lying non-zero modes of the overlap Dirac operator which 
seem to accrue to topological structures and exhibit strong chirality, as measured by the local parameter X defined 
by 



tan 



in regions where the probability density ijj^ (x)tp(x) of these modes is maximal ||. The verification of the instantonic 
nature of these objects and their relevance to spontaneous chiral symmetry breaking in the infinite volume is still 
being argued out in the literature (see for example the recent large N c study of ||). 

An unbiased summary of the totality of available lattice results can be formulated as follows: that they support the 
importance of gluon configurations producing regions of approximately "locked" chromo-electric/magnetic fields for 
chiral symmetry breaking but do not yet confirm or rule out a specifically instantonic nature for these configurations |j) . 
A potential test which might clarify this would be a comparison of hadronic correlation functions between vacuum 
models and lattice simulations. Such results are already available for instanton-based models JhJ. A search for 
complementary scenarios for the vacuum consistent with lattice results and incorporating confinement (missed in the 
instanton models) is evidently timely. This paper is a step in that direction for the domain-scenario. 

The core of this paper is the Dirac eigenvalue/function problem for a spherical four-dimensional Euclidean region 
of radius R with bag-like boundary conditions on the fermions and in the presence of a covariantly constant (anti-) 
self-dual gauge field 

P1>(x) = \i/>(x), 

i fi{x)e ia ^i){x) = i/>(x), x 2 = R 2 . (2) 
Here rj^(x) = Xfi/\x\, is the covariant derivative in the fundamental representation, 

= 9 M - iB^ = d^ + -hB^Xv, 
h = t a n a , = ±B IW , 

where the (anti-)self-dual tensor £? M „ is constant, and the Euclidean 7— matrices are in an anti-hermitean representa- 
tion. 

The outcomes of this study are the peculiar chiral properties of eigenspinors ip(x): there are no zero modes and 
none of the modes is chiral but at the centre of a domain the local chirality parameter X{x) is found to be 

X(0) = ±1 

for all modes with zero orbital momentum. The sign of chirality and the duality of the tensor B^ v are locked: 
(+1) — 1 for an (anti-)self-dual field. Simultaneously the normal density for these modes is maximal at the centre. 
At the boundary the local chirality X is equal to zero for all modes. The chirality of the lowest mode is a monotonic 
function inside the region while for the higher radial excitations the chirality alternates. The detailed form of X(x) 
changes with the variation of an arbitrary angle a in the boundary condition. This angle is treated as a random 
variable. Calculating chiralities averaged over a small central region for the various lowest modes, and operating 
in the whole ensemble of domains, we end up with a histogram which represents the probability of finding a given 
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smeared chirality among the set of lowest modes. The histogram qualitatively reproduces the lattice results for the 
chirality of low-lying Dirac modes such as those of || and others. 

After reviewing the domain model in the next section, we present details of the solution of the above-formulated 
problem and then study the chirality properties of the eigcnmodes. We conclude with a discussion and future prospects. 
Technical details of calculations and conventions for this paper are relegated to the Appendices. 

II. REVIEW OF THE MODEL 

It has been suggested (TTJ that the restrictive influence of pure gauge singularities (present in instanton, monopole 
and vortex configurations) on surrounding quantum fluctuations may be used for an approximate treatment of QCD 
dynamics. Due to the complex structure of the manifold of gauge orbits in QCD, singular gauge fields may be 
unavoidable in the course of the elimination of redundant variables. Obstructions such as the Gribov problem and 
condensation of monopoles are two examples of this potentially more general statement. This has also long been 
advocated by van Baal in his studies of the fundamental domain in small volume studies on the torus and sphere. 
In particular, the proposal has been made that "domain formation" at larger volumes can be driven essentially by 
the non-trivial topology of the gauge field manifold. Moreover, it is stressed in Q that the full set of singular 
fields, instantons, monopoles and vortices must play a role in this. One can also add to this hierarchy domain wall 
singularities [ p"3| which are not topologically stable on their own but can be part of a complicated object: a domain 
wall can start and end on a lower dimensional topologically nontrivial singularity of lower dimension, namely a vortex, 
and in this sense should not be neglected also. 

An arbitrary gauge field configuration A containing a pure gauge singularity can be represented in the vicinity of 
the singularity as 

•A^i — -t- Q [_l 

with Sfj, a pure gauge singular field. If we now substitute this into the Yang-Mills Lagrangian we will see that the 
requirement of finitness of the action density imposes specific conditions on the behaviour of Q in the vicinity of the 
singularity in S. The model we consider focusses on domain wall singular hypersurfaces which are the most restrictive 
for Q; an inclusion of lower dimensional singularities is a complicated task beyond the scope of the present work. In 
the case of domain wall the constraining influence of S on gluon fluctuations Q and quark fields ip is expressed via 
the boundary conditions 

[Q,S} = 0, (3) 
t/j(x) rj(x)ip(x) = (4) 

for x being on the singular hypersurface of the pure gauge field S. These conditions ensure a non-vanishing weight 
for such fields in the functional integral. 

Domain wall singular pure gauge configurations are topologically trivial. This implies that the field S can be 
characterised by a definite colour direction n a and the matrix n a t a can always be tuned to belong to the Cartan 
subalgebra of SU C (2>). The off-diagonal (or, equivalently, orthogonal to n a ) components of the fluctuations Q must 
then satisfy Dirichlct boundary conditions, while those fluctuations longitudinal to n a arc not restricted at the domain 
wall. A typical configuration of this type looks like a system of domains which are coupled in a sense that fluctuations 
inside neighbouring domains interact with each other via exchange by the gluon modes longitudinal to the colour 
direction of the domain boundaries. It should be stressed that unavoidably there are obstructions of colour direction 
at the domain wall junctions where lower dimensional topologically nontrivial singularities are situated. 

To be specific and to deal with an analytically tractable model we introduce several drastic simplifications: we 
disengage ourselves from the obstructions in the colour direction and substitute the coupling between domains by 
the presence of a mean field. Inside and on the boundary of the domain the field is taken to be covariantly constant 
(anti-)self-dual such that the strength over the whole Euclidean space reads 

N 

i=i 

The individual colour and space orientations in each domain are random. In particular, effective action arguments 
were used in |l| to constrain the form of rS^ a such that the matrix = i 3 cos£j + i 8 sin£j with angles £j S 
{■f (2fc + 1), k = 0, . . . , 5} corresponding to the discrete Weyl subgroup. Domains are taken to be hyperspherical with 
a mean radius R and centered at random points Zj . For a detailed motivation of these steps we refer the reader to [fy . 
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In this way, consideration is reduced to a model with essentially two free parameters: the mean field strength B 
and the mean domain radius R. The partition function for this simplified system can be written down as 

N 

Z = M hm T] daj ThfiWvj® / 2>Q\J[D(tf«)QM] AFp [g« Q W ]e -s$° D [Q w +B™ ,1>™,i> w ] (6) 

where the functional spaces of integration Tq and T\ are specified by the boundary conditions (x — Zi) 2 = R 2 

h l Q^(x)=0, (7) 
ifn{x)e ia ™i)^{x)=i)^{x), (8) 
^e^ifj^x) = -$M(x). (9) 

Here fii = nft a with the colour generators t a in the adjoint representation. The conditions Eqs.(^) and (||) represent 
specific (though not unique) choices for the implementation of Eq. (|J) which manifest the explicit breaking of chiral 
symmetry by the boundary condition, as occurs for example in bag models for the nucleon. The thermodynamic limit 
assumes V, N — > oo but with the density v -1 = N/V taken fixed and finite. The partition function is formulated 
in a background field gauge with respect to the domain mean field. The measure of integration over parameters 
characterising domains is 



/ 



2tt 

d(J i--- = T3~T / —rr /) dcti \ difi / dOi sin 0, 

487r 2 ./ y V ^^Jpvi-i)* J Jo 



E — * 1 

x r ^ e - (2 ^t i)7r ) r e - ?rfe ) ■ ■ ■ . ( iq ) 

J ° 1=0X2 ° J ° fe=0,l 

where (6i,tpi) are the spherical angles of the chromomagnetic field, u>i is the angle between chromoelectric and 
chromomagnctic fields and is an angle parametrising the colour orientation. It should be noted that because of 
the axial anomaly and that nothing a priori constrains the topological charge per domain to be integral the fermion 
determinant is a single valued function of at only if an appropriate Ricmman surface is constructed. Here i>i enumerates 
the Riemman sheets to be taken into account. 

This partition function describes a statistical system of density v~ l composed of extended domain-like structures, 
each of which is characterised by a set of internal parameters and whose internal dynamics are represented by the 
fluctuation fields. It respects all the symmetries of the QCD Lagrangian, since the statistical ensemble is invariant 
under space-time and colour gauge symmetries. For the same reason, if the quarks are massless then the chiral 
invariance is respected. 

Field eigenmodes satisfying the above boundary conditions in the presence of an (anti-)self-dual gluon field and 
corresponding Green functions can be calculated explicitly. For gluons, this was shown in [jjj. For quarks, this will 
be shown in this paper. On this basis one can compute any correlation function taking the mean field into account 
exactly and decomposing the integrand over fluctuations. In particular, correlation functions of the mean field itself 
have a finite radius i?, which is more or less obvious and is discussed in |l) in detail. 

Within this framework the gluon condensate to lowest order in fluctuations is immediately obtained in the form 

9 2 {f; v {x)f; v {x))=^b 2 , (ii) 

and the topological susceptibility reads 

B 4 R A 



X = / d i x(Q(x)Q(0)) 



1287T 2 



Less trivial is the manifestation of an area law for static quarks. Computation of the Wilson loop for a circular contour 
of a large radius L^> R gives a string tension a = B f(TrBR 2 ) with the function 

2 / ^3 r 2z/v * 




Estimations of the values of these quantities are known from lattice calculation or phenomenological approaches and 
can be used to fit B and R. As described in Q these parameters are fixed to be 

VB = 947MeV, R = (760MeV) _1 = 0.26fm (12) 
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with the average absolute value of topological charge per domain turning out to be q « 0.15 and the density of 
domains v^ 1 = 42fm~ 4 . The topological susceptibility then turns out to be x ~ (197MeV) 4 , comparable to the 
Witten-Veneziano value Q. 



III. SPECTRUM OF THE DIRAC OPERATOR IN A DOMAIN 



We have mentioned already that the boundary conditions on fermions violate chiral symmetry explicitly which can 
only be restored by a random assignment of values of a over the complete ensemble of domains in Euclidean space. 

In this section we address the eigenvalue problem for the massless Dirac operator as it is stated in Eqs. (^). Here 
we give the scheme for solving the problem, with technical details given in the appendices. The Dirac matrices in 
Euclidean space are chosen to be anti-hermitean and taken in the chiral representation. 

For boundary conditions on a hypersphere and covariantly constant background field of definite duality it is natural 
to use hyperspherical coordinates (r, f2), given in detail in Appendix A. In such coordinates, rather than work with 
the covariant derivative itself, it is more convenient to introduce the operator fj p which can be easily expanded into 
intrinsic and orbital angular momentum generators. Any spinor can be represented in the form 

ip = iyjx + <p, ip = ixyj+(p, (13) 

where tp and x have the same chirality. This is simply a decomposition into a sum of chiral components. The 
eigenvalue equation Eq. (|^) can be rewritten then identically as 

x = ~fiPf, P\ = \ 2 y. (14) 

In these terms the boundary conditions take the form 

X = -e T "V, x = Ve Tia , x 2 = R\ (15) 
where upper (lower) signs correspond to (p and x with chirality =pl. 



A solution of Eqs. (14) is achieved by separating the angular and radial coordinates. To do this one has to represent 
respectively p 2 and fjp in terms of momentum generators and projectors onto the various spin and colour polarization 
subspaces. In four-dimensional Euclidean space the angular momentum operators can be represented as 

Ki, 2 = i(L±M) 

with L the usual three-dimensional angular momentum operator and M the Euclidean version of the boost operator. 
These correspond to the decomposition of the four dimensional rotational group SO(A) into a product of two 5*0(3) 
groups. They lead to Casimirs and eigenvalues 

K? = K»-|(£ + l), fc = 0,l,...,co 
#1,2 -> mi, 2 , TOi, 2 = -k/2, -k/2 + 1,... ,k/2 - 1, k/2, 

and the correponding angular eigenfunctions Cfe mim2 (0), given explicitly in Appendix A, are labelled by orbital mo- 
mentum k and two azimuthal numbers mi and 7712. Eigenstates are also characterized by the colour-spin polarisation 
related to the projectors 



0±=A + S ± +A_E T (16) 



with 



N±= 1 -{l±n/\h\), E± = i(l±SB/B) 

being respectively the separate projectors for colour and spin polarizations. Below we denote the polarisation with 
respect to O by K = ±. 

It is shown in Appendix B that if the background field is (anti-)self-dual the boundary condition can only be 
implemented if spinors (p and x are (right) left handed. Also the presence of the homogeneous background field 
reduces the spherical symmetry of the problem down to an axial-symmetry. In the representation implemented here 
this manifests itself as a restriction on the values of one of the azimuthal quantum numbers, namely to 2 = ±fc/2 
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for the self-dual case and mi = ±fc/2 for the anti-self-dual one. The sign in front of fc/2 is correlated with the spin 
polarization of the state as seen in the explicit expressions for the eigenspinors below. 

Thus for the self-dual case, j^ip = —ip, 75X = — X so that the eigenspinors in the self-dual field can be labelled as 
ipkm {x), while in the anti-self-dual field they are i'Xm.J^)- With details in Appendix B, we simply write down here 
the result for the self-dual case: 



' femi 



-(iA)- 1 z( fe+1 )/ 2 e-^ 2 [M (k + 2 - A 2 , k + 2, z) 

( 

k + 2 - A 2 



(17) 



fc + 2 



-M(fc + 3- A 2 ,fc + 3,z) 







^ = z k / 2 e- z / 2 M(k + 2-k\k + 2,z) 



\N + c kmi _,(n) J 

° \ 



z (k+l)/2 e -z/2 



iA 
fc + 2 



M (l - A 2 ,fc + 3, 



/ 



? V2 e -*/2M(-A 2 ,fc + 2,. 



N-C kmi ,(n) 
\ N+c kmi _,(n) J 

I 



N+C kmii (Q) 

2 ^ 



(18) 



(19) 



where M(a,b,x) is the confluent hypergeometric function and 

z = Br 2 /2, A = X/V2S, B = \h\B. 

The projectors N± act on the colour vectors which are implicit in above equations. The eigenf unctions V , fc~m 2 f° r the 
anti-self-dual case are obtained by the change mi — > m2 and the shift of nonzero elements of the angular part to the 
first two positions of the spinor. 

The eigenvalues are determined by the boundary condition at z = zq = BR 2 /2 which for A^ + takes the form 

k + 2 - A 2 



e~' la M (k + 2 - A 2 , k + 2, Zq) - ^ 

iA 



M (k + 2- A 2 ,fc + 2,z 



-M(k + 3- A 2 ,fc + 3,z ) 



and for A, 



e - m M(-A 2 ,fc + 2,2 ) 



iA 



M(l 



A 2 ,fc + 3,z ) = 0. 



0, (20) 



(21) 



The equations for the eigenvalues in an anti-self-dual domain are the same as above but with a — > —a as follows 
from Eqs. (p"5|). The eigenvalues can be calculated numerically. They form a discrete set. Zero modes are absent, 
which is to be expected for these types of boundary conditions fl5|| . A graphical solution of Eq. ( ^l|) at a = ir/2 is 
presented in Fig.|l| to illustrate the structure of the spectrum. In general, the eigenvalues are complex. The spectrum 
is real for a = ±ir/2 which is the only value for which the boundary condition Eq.(Q) imposed on ip^^i 2 * s nerm itean 
conjugated to the condition for il>km 2 an< ^ ^ e § enera l fermion field ip can be decomposed in terms of the basis of 
conjugate eigenfunctions 4> k ^ ll ■ For other values of a a biorthogonal basis should be introduced. In particular at 
a = eigenvalues are complex and come in complex conjugated pairs. The partition function is nevertheless real 
since if A s d(aO is an eigenvalue for the self-dual case then for the anti-self-dual domain there is an eigenvalue A as d(a) 
such that 



Aasd(a) = -A* d (a). 



(22) 



We stress that the definition of A here does not include an imaginary unity in front of Dirac operator. 
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As seen from Fig.[I], in contradistinction to the eigenvalue problem in infinite volume on the space of square integrable 
functions the spectrum is not symmetric under reflections A — > —A. This comes from the fact that 75 does not 
commute with the boundary condition so that 75^ is not an eigenfunction if ip is an eigenfunction. An assymmetry 
of the spectrum is typical for the Dirac operator in odd dimensional spaces (see |]l6f and references therein) and has 
important consequences there for the effective action. In our case the unusual boundary conditions are responsible 
for the assymmetry in four dimensional Euclidean space [ fL5| . 

The most interesting feature of the fermionic eigenmodes becomes manifest if one considers the local chirality X(x) 
of the lowest eigenmodes as defined by Eq. (|T]). 




FIG. 2: Normal density distribution for modes with k — and n — 1, 2, 3. 



IV. CHIRALITY OF LOW LYING MODES 



It is obvious that none of the solutions are eigenstates of 75. However at the domain center x 2 = (or (x — Zj) 2 = 
in general) all the purely radial modes with k = have a maximum in the probability density, they are chiral and 
the sign of their chirality is determined by the duality of the mean field in a domain which is illustrated in Figs. 2 and 
0. The probability density naturally vanishes at the domain center for the modes with k > as is seen in Fig. 4 . To 
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Distance from the center (fm) 



FIG. 3: Chirality parameter for the three lowest radial modes %p 00 , self-dual domain, a — tt/2. 



demonstrate this analytically let us turn to the local chirality parameter given in the introduction which we rewrite 
here in more detailed form 



tan x( X ))) = - 75 ^ , -i < x ix) < i, 

which takes the extremal values X = ±1 at positions x where ip(x) is purely right(left) handed. Because ip and x 
have the same chirality, the representation Eq.(|l7|) immediately gives for the self-dual domain 

taa(J(l-X-/(,)))=l^M, 
V4 ' IXoo 0)1 

while for the anti-self-dual case the local chirality reads 

ta n (J(l-X-(,)))=^gM. 
V4 ' ma \ x )\ 

Moreover due to the relation Eq.(|22]) 

l<4f(*)l = W(*)l> IxM = IXoTOI- 
Representations Eqs. ([lj) and ( |l9| ) show that 

< lira lv>ou»l < oo, Km \x$o(x)\ = 0, 

which finally results in 

x+/(o) = -i, x+2(o) = i. 

The local chirality parameter X as a function of distance from the domain centre for the lowest few modes is plotted 
m Fig§ There is a peak in X at the domain centre. Away from the centre X decreases due to a competition of 
left and right components of the eigenmodes as the x component becomes non-vanishing. As is seen from Fig.|| the 
chirality of the lowest mode (n = 1) monotonically decreases with distance from the centre. The chirality parameter 
for the excited modes alternates between extremal values, the number of alternations is correlated with the radial 
number n, and the half-width decreases with growing n. The chirality parameter X is zero at the boundary for all 
modes. Qualitatively this picture does not depend on the angle a. The "width" of the peaks at half-maximum for 
the lowest (n — 0) radial modes varies for different values a and is of the order of 0.12 — 0.14 fm if the values of B 
and R are fixed from the gluon condensate and the string tension, consistent with the lattice observations of Q. 

We can now study the chirality characteristics of the ensemble of fermion fields entering the partition function 
Eq. ( 3) with all values of a treated with equal probability consistent with an explicit chiral symmetry. On the lattice 
H |?] peaks in X or ip^ip would only be localiseable within a size corresponding to the lattice spacing. To take this 
into account, we average X(x) over a small neighbourhood of the domain centre. Thus we compute the probability 
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to find a given value of X, the smeared X, among the chiralities for the lowest modes. The result given in Fig.|5| 
was obtained for three sets of modes: with n < 2, k = (solid line), n < 4, k — (dashed line) and n < 6, k = 
(dot-dashed line), and all possible values of a and spin-colour polarizations. The solid line, formed from the lowest 
modes, evidently indicates two narrow peaks with X ±0.87. This double peaking is not unexpected in view of 
the above discussed chirality properties. Including higher modes broadens the peaks and shifts their maxima. This 
feature as well as the above mentioned values for the half-width and the density of domains is in qualitative and 
quantitative agreement with recent lattice results ||, |l7|]. It should be stressed that orbital excitations (k > 0) are 
not included in the histograms, because the probability density for orbital modes vanishes at the centre. However 
there are maxima in the probability density for these modes in peripheral regions of the domain. The local chirality 
X is significantly smaller in peak value than those for the radial modes at the centre. Inclusion of orbital modes will 
broaden the peaks more and build up the central plateau. 




FIG. 4: Plot of the radial dependence of the normal density distribution for the modes with k — 1 and n = 1,3, 5. 
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FIG. 5: Histogram of chirality parameter X averaged over the central region with radius 0.025fm. Plots given in solid, dashed 
and dot-dashed lines incorporate all modes with n < 2, n < 4 and n < 6 respectively. 
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V. DISCUSSION AND CONCLUSIONS 

The statement that signals for spontaneous chiral symmetry breaking should be identifiable in the specific chirality 
properties of fermionic eigenmodes for some "dominant" gluonic background field is generally accepted. Such signals 
are now being seen on the lattice, but nevertheless there are not so many analytically explicit examples of this 
relationship available. Instanton motivated models are certainly the most advanced example of this kind. But one 
aspect of instanton models lies in the 'dual' use to which the exact fermion zero modes of the classical instanton 
background are subject: these modes have definite chirality and do not contribute in the infinite volume limit, but 
are used to model pseudo zero modes which should characterise the non-semi-classical system of strongly overlapping 
instantons and anti-instantons. 

We have studied the spectrum of quark modes in a domainlike structured gluon background field. Such a background 
is argued to characterise the bulk average properties of the vacuum in the presence of strong intermediate range 
fluctuations and is not the result of a semi-classical approximation. The spectrum exhibits definite chirality properties. 
In particular, there are no zero modes because of the conditions which fermion fields must satisfy on the boundaries 
of domains. Nonetheless at the centre of domains all radial modes are purely chiral and the sign of their chirality 
depends on whether the underlying gluon field is self-dual or anti-self-dual. Moreover, the sign of chirality at the 
centre persists over the whole domain for the lowest modes. Studying the local chirality parameter X in a chirally 
symmetric ensemble of domains we obtain qualitatively similar results to those seen in lattice calculations. We stress 
that this comparison with lattice results takes place at the level of an ensemble of configurations not on a configuration 
by configuration basis. 

Insofar as these lattice results for chirality are argued as supporting the evidence for spontaneous chiral symmetry 
breaking, the same can be said of the domain model. We note the absence of any explicit zero modes in achieving this. 
Namely, the range of configurations needed to produce the types of effects seen on the lattice are not restricted to 
instanton-like fields. It suffices that a given gluon background admit strongly chiral low-lying non-zero modes. In this 
respect, the more significant property of the gluon background is the "locking" of chromoelectric and chromomagnetic 
fields into self-dual or anti-self-dual fields in relatively large but finite regions of space restricted by the hypersurfaces 
on which a pure gauge singularities are assumed to be situated. It should be stressed that in the thermodynamic limit 
the number of domains is growing but their sizes stay fixed around some finite mean value. 

The solutions obtained in this paper provide a basis for computation of chiral condensate (ipijj) , in particular in the 
presence of an explicit CP violating 9 term. In view of the assymmetry of the spectrum, the specific chiral properties 
of the eigenmodes and noninteger mean topological charge associated with a domain are expected to simultaneously 
enhance the role of the axial anomaly in the spontaneous breakdown of chiral flavour symmetry and factor in the 
solution of the I7 A (1) problem || [[§, |§. 
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APPENDIX A: CONVENTIONS 

We use a chiral representation for the antihermitian Dirac matrices in four Euclidean space, 

{it,, lA = -2<5 M „, 7^ = 

7i= (-cr, )' 74 = Z (5 o)' 
75 = 71727374 = diag(l, 1, -1, -1). 

The background field is specified as 

Bij = EijkBk, Bi — -SijkBjk, Ei = Bi4 = ±B i: 
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In addition the following conventions and relations have been used 

1 _ 1 

CijBij — 2T,iBi, E± = — f 1 ± 



and in particular 



B 

1 1 



Va&B a p = (JijBij + 2aaBa = 2S^i?i =F 2J5i7sSj 
= 4P T E i B i = 4BP q :(E + -S_). 



We use the following hyperspherical coordinate system in R 4 



Xi = r sin 77 cos 0, x 2 — r sin 77 sin 0, 
x 3 = r cos 77 cos x, £4 = r cos 77 sin x, 



(Al) 



and define the angular momentum operators as follows 



Mi = -i(x 4 di - Xidi), 

respectively for spatial rotations and Euclidean "boosts" . As mentioned, it is more convenient to work in the basis 

Ki >2 = (L±M)/2 

which generates the following Lie algebra 



K\,K{ 



KlKi 



= ie ijk Kk, 



Thus the ladder operators 
satisfy the algebra 



= ie ijk Kj[, 

[Kl 2 ,K± 2 ]=±K± 2 



K[,Kl 



= 0. 



and correspond to raising and lowering operators of mi , m 2 . 

The angular eigenfunctions corresponding to the Ki j2 generators are 

C kmim2 (v,<P,x) = (-1)1^+^1 (27r)- 1 er~ m2 ' mi+m2 W, 
x cxpi[(mi — m 2 )x + (mi + m 2 )4>] 

e k k - r - s ' s - r (r]) = y/2(k + l)(k - r)\(k - s)\r\s\ 

' (-l) r - n cos k - r - s+2n ri sin r+s - 2n rj 



(h — r — 



£-^ Q (k — r — s + n)\n\(r — n)\(s — n)\ 
s = (k + m 2 )/2, r = (k- mi)/2, 

where k,nii,m 2 are respectively the orbital angular momentum and the two azimuthal quantum numbers, relevant 
for an 0(4) = (9(3) x (9(3) symmetry. They take the following values 



k = 0, 1,2, . . ., mi, 7712 



k k 

' o ' ' ' * ' n " 
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APPENDIX B: DIRAC EIGENVALUE PROBLEM IN A DOMAIN. 



Here we give further details of the solution of Eq. (Q) . Using the notation given in the main body, we can decompose 
the field if over a set of chiral and colour-spin projectors 



P±$ 



P t O+$+i 



where (lower)upper signs correspond to the (anti-)self-dual field background field, and fields <£>£ must satisfy the 
second order equation 



-D 2 + 2CB 



0. 



(Bl) 



We remind that implicitly <!>" is the color vector in the fundamental representation. 

If we were solving the problem for square integrable eigenfunctions in in finit e volume then all three components 
would enter the final set of eigenfunctions, moreover as is seen from Eq.(Bl), the equation for the component $_i 
would produce zero modes with chirality =pl. The spectrum in this case would be discrete and all nonzero eigenvalues 
of Dirac operator come in pairs: if ip is an eigen function with eigenvalue A then 75-0 is an eigenfunction with eigen 
value —A - pretty standard state of affairs. 

However the bag-like boundary conditions Eq.(|^) we must satisfy in the present case change the structure of 
eigenfunctions and eigenvalues drastically. First of all, because of the identities 



"j^B^pXpP^Y,. 



= iB ^P T S+, 
= -iB ^P T S_ 



and 



jpB t _ t px p p ± j: + 

r yp,Bp P XpP±Y,_ 



2i 



B[x 2 - [B^f/B 2 ] - iB \ ^P±£+ 



2 1 



B[x 2 - {B lXl ) 2 /B 2 } - iB \ ^P±E 



which can be straightforwardly derived by expanding both sides over a complete set of Dirac matrices, the bag-like 
boundary condition can be satisfied only for the trivial solution &o(x) = 0. The significance of this observation is 
that for (anti-)self-dual domains the boundary condition can only be implemented on eigenspinors tp = i fix + 4> for 
(positive) negative chirality if and x- The function ip in turn is not an eigenspinor of 75 , which is natural because the 
boundary condition violates chiral symmetry. Furthermore, zero modes are removed from the spectrum because they 
must be chiral but this is forbidden by boundary conditions. And, finally, if ip is an eigenfunction with eigenvalue A 
then 75^ is not an eigenfunction any more, and there is no eigenvalue —A in the spectrum. 
In order to find equations for components of the corresponding spinors we use that 



f-± =P-0±*±i - (0,0,7V T $ 3 ±1 ,iV ± <i> 4 ±1 ) r 
f+± =P+0±*±i = (A f T $i 1 ,iV ± $| 1 ,0,0) r 



(B2) 



In hyperspherical coordinates Eqs.(Al) the equations for the spinor components read (here and below we write down 
equations for the self-dual case only) 



7 d r r 3 d r 



rd r r 3 d r 



—d r r 3 d r 



-K 2 



Tl n 

2—BK 2z 
\n\ 

h - 
2-~BK 2z 
\n\ 

71 * 

2--BK 2z 
\n\ 

h ~ 
2—BK 2z 
n 



-B 2 r 2 



B 2 r 2 



-B 2 r 2 



-B 2 r 2 



2B - X > iV_$ 



The anti-self-dual case is reconstructed by the change K 2z — > K\ z and $ 



C 



2B - X 



S3 



0, 



N+^ +l = 0, 



2P - A > JV+$ 



0, 



— 2B — X > N-§\ = 0, 



$1, $4 



$ 2 . 



general solution for equations of this type. The requirement of regularity at the origin then gives 



(B3) 



In |l| we derived the 



^ mim2 = N„z k / 2 e- z / 2 M(t+m 2 -A 2 + 2,k + 2,z)c kmim2 ( V >,X,v), 
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$^ mim2 = iV + z fe/2 e - z/2 M^-m 2 -A 2 + 2 ; fc + 2,z^C femim2 (^ X ,7 7 ), 
= N + z k ' 2 e- z l 2 M Q - m 2 - A 2 , k + 2, ^ C fcmi „ l2 fo>, X , »?), 
= N_z k ' 2 e- Z I 2 M (^+m 2 -A 2 ,k + 2,z^jC kmim2 (^x,v), 



( j ) 4,fcmirri2 



where A = A/ \/ 2B, z = Br 2 /2. Thus the two independent mutually orthogonal solutions 



arc 



-+ _ 








\ 



N _ z k/2 e -z/2 M (| + m2 _ A 2 + 2,fc + 2, z) C femim2 
v 7V + z fe '/2 e -,/2 M (V _ m / 2 _ A /2 + 2 , fc' + 2, z) C fe , m ; m / y 



<y5 = 



\ 



N + z k/2 e~ z/2 M (| - m 2 - A 2 , A; + 2, z) C 



krn\m,2 



(B4) 



(B5) 



where "prime" indicates that angular quantum numbers and eigenvalues in the third line need not coincide with those 
in the fourth in order that these spinors be eigenmodes of Eq.(|l4|). 
To obtain an explicit representation for \ we use the identity: 



- ff(x) p = d r + 2R-\-E • KxP + + S ■ K 2 P_) - i ff^B^x v , 



(B6) 



where the action of the last term on PzpO^Q** can be determined via the identity 



fi BR 
-l^B^XvP^Oc, = i(fj—P^O c , 

and the action of the £ • K terms via 

e.k 1)2 £ o c $f = (E 3 xr i2 + s (+) ^r: 2 + s(-)x+ 2 ))^o c $ c 

c c 

|n] 

+iV + E( + )if 1 - 2 $ +1 + N-& +) Ki2$-i + N_Y±(-)k£ 2 $ +1 + N + J: { - ) K+ 2 <P^ 1 . (B7) 
As well as the ladder operators K ± we also have analogous operators for the spin, 

E (±) = £(Ei±i£ 2 ). 



For i?i = -6(5^3 the following identities are also useful for implementing the above 

n 



and 



P_E (+) * = (0,0,-* 4 ,0) T , 

p_s(-)* = (o,o,o,-* 3 ) T , 

P+£ (+ )* = (-* 2 ,0,0,0) T , 

P + E (_) tf = (0,-W 1 ,0,0) T . 



(B8) 
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We thus get for the self-dual case 
/ 



and 



j-.[d r -Br_ l K l\N-z k l 2 e- z / 2 M (| + to 2 - A 2 
Mdr -^f + iKl]N + z k '/ 2 e-*/ 2 M 



( 



A 



2, -|- 2, Cf i - mim2 

- 2, -f- 2, J Ck'ra', mL 



\ 



■\jKzN + z k '/ 2 e- z l 2 M ( Hr ~ to' - A' 2 + 2, jfc' + 2, z ) C fe , 



-'k'm 1 -. mi. 



\ ^K+N_z k / 2 e-^ 2 M(l+m 2 -A 2 + 2,k + 2,z)C kmim2 J 



( 



\ 



Txl d r + ¥~ $Kl]N + z k / 2 e- z / 2 M (f - m 2 - A 2 , k + 2, z) C kmim2 
{ M d r + ¥ + %Ki]N_z k '/ 2 e-*/ 2 M + m' 2 - A' 2 ,k' + 2, z) C fe , TO , m , y 

/ 



I i 7 X 2 -iV_z fe '/ 2 e- z / 2 Af (| + to 2 - A' 2 , fc' + 2, zj C fc , m; „ 
^ ±K+N + z k l 2 e- z ' 2 M (| - m 2 - A 2 , + 2, z) C femim2 



\ 



(B9) 



(BIO) 



By inspection, the boundary condition \ — ~e~ la ip can only be fulfilled if terms with raising/lowering operators of the 
azimuthal quantum numbers vanish since t hes e term s co ntain the projectors N± while the rest of terms entering the 
boundary condition contain N T (see Eqs. (B4) and (JB5|)). In particular, m' 2 — —k/2, to 2 = k/2. Finally, evaluating 
the derivatives of the confluent hypergeometric functions with the help of relation 

M' g (a,b,z) = ^M{a + l,b+l,z) 

leads to the solutions given in the main body of the paper. 



[1] 
[2] 



[3] 

[4] 

[5] 
[6] 
['] 
[8] 
[9] 
[10] 

[11] 

[12] 
[13] 
[14] 
[15] 
[16] 
[17] 
[18] 

[19] 



A.C. Kalloniatis, S. N. Nedelko, Phys.Rev. D 64, 114025 (2001). 

Some selected references are: M.J. Teper, Phys. Lett. B162, 357 (1985); E.-M. Ilgenfritz et al., Nucl. Phys. B268, 693 
(1986); M.J. Teper, Nucl. Phys. B (Proc. Suppl.) 83, 146 (2000); UKQCD Collaboration, Douglas A. S mith et al, Phys 
Rev. D 58, 014505 (1998); Philippe de Forcrand et al., Nucl. Phys. B (Proc. Suppl.) 63, 549 (1998) see als o bcp-lat/9802017 . 
For a relationship between cooling and smearing see F. Bonnet et al., Phys. Rev. D 62, 094509 (2000); hep-lat/0106023. 
T.L. Ivanenko and J.W. Negele, Nucl. Phys. B (Proc. Suppl.) 63, 504 (1998); R.G. Edwards, U.M. Heller, and R. Narayanan, 
Nucl. Phys. B522, 285 (1998). 

J.B. Zhang, S.O. Bilson-Thompson, F.D.R. Bonnet, D.B. Leinweber,A.G. Williams, J.M. Zanotti, Phys.Rev.D 65, 074510 
(2002). 

Robert G. Edwards, Urs M. Heller, Phys.Rev. D 65, 014505 (2002). 
T. Banks, A. Casher, Nucl.Phys.B167, 215 (1980). 

I.Horvath, N. Isgur, J.McCune, H B. Thacker, Phy s.Rev. D 65, 014502 (2002). 
N. Cundy, M. Teper, U. Wenger, |hepdat/0203030 



R.G. Edwards, Nucl.Phys. B (Proc.Suppl.) 106, 38 (2002). 
T. DeGrand, Phys.Rev. D 64, 094508 (2001). 

O. Jahn, F. Lenz, Phys. Rev. D 58, 085006 (1998); F. Lenz S. Worlen, At the frontier of particle physics, vol. 2, 683-760 
ed. by M. Shifman, World Scientific, 2001 ( |hep-ph/0010099[ ). 
P. van Baal, At the frontier of particle physics, vol. 2, 760-821, ed. by M 
C. Ford et al, Ann. Phys. (N.Y.) 269, 26 (1998). 



P. Di Vecchia, G. Veneziano, Nucl.Phys.B171 
A. Wipf and S. Duerr, Nucl. Phys. B443, 201 
S. Deser et al, P hys. Rev. D 57 7444 (1998). 
I. Horvath et al. bep-lat/020100£ . 
R.J. Crewther, "Chiral Properties of 
(QCD161:N16:1979); Nucl.Phys. B209, 413 (1982) 
R. Dashen, Phys. Rev. D 3, 1879 (1971). 



253 (1980); 
(1985). 



E. 



Shifman, World Scientific, 2001 ( |hep-ph/0008206[ ) 
Witten, Annals Phys.128, 363 (1980). 



Quantum Chromodynamics" , Kaiserslautern NATO Inst. 1979:0529 
See also, R.J. Crewther, Phys.Lett.B70, 349 (1977). 



